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Abstract 

One of the most powerful tools to probe the existence of cosmic defects in the early universe is 
through the Cosmic Microwave Background (CMB) radiation. It is well known that computa- 
tions with causal sources are more involved than the adiabatic counterparts based on inflation, 
and this fact has in part hampered the development of fine detailed predictions. Analytical 
modeling, while necessarily limited in power, may tell us the overall characteristics of CMB 
from defects and hint at new features. We apply an analytical model for textures to the study 
of non-Gaussian features of the CMB sky and compare our predictions with the four-year 
CttB^DMR data. 
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1 Introduction 



It has by now become clear that one of the most promising ways to learn about the early universe 
is through the Cosmic Microwave Background (CMB) radiation. With the prospective launch 
of future missions, like MAP and COBRAS/SAMBAQ, one can hope that many of the so far 
elusive cosmological parameters will be pinned down with unprecedent precision. The four- 
year COBE plus other large-scale structure data have placed strong constraints on current 
models of structure formation. However the remaining window is still too large and many 
(widely) different cosmological models pass the test. This is actually the case with cosmic 
defect theories, on the one hand and inflation-based adiabatic models on the other. 

The search for the so-called Doppler peaks (or Sakharov peaks if you like) seems to be 
among the first goals for next generation detectors, the aim being in trying to discriminate 
among, say, standard adiabatic CDM models 1 ^, cosmic strings 2 -* and textures 3 - 1 . Whereas the 
generation of CMB anisotropies seems to be fully understood within adiabatic models [refer to 
Hu's contribution to these proceedings] , the same does not happen for the latter models, where 
the non-linear evolution of the defects and their active role in seeding anisotropies in the CMB 
makes the analysis far more involved [see Durrer's contribution]. Moreover, causal sources 
can produce spectra mimicking the outcome of inflationary models 4 -*, hence increasing the 
uncertainty and calling for fast and accurate methods for computing the theoretical predictions 
[Seljak's contribution], and refined experiments to confront these competing theories. 

Other means of narrowing somewhat the window regards the recognition that the CMB 
radiation carries valuable information of the processes that generated the anisotropies, all along 
the path of the photons from the last scattering surface to our present detectors: should future 
experiments find (with high confidence level) a departure of the statistics of the anisotropies 
from Gaussianity, it would disfavor the standard inflaton field quantum fluctuation origin of the 
cosmological perturbations. It then follows that it is interesting to calculate what predictions 
cosmic defects make regarding non-Gaussian features in the CMB sky. It is the aim of this 
short contribution to report on some work done on the CMB three-point correlation function 
predicted by textures within a simple analytical model. 



2 The model 

Magueijo 5 ) recently proposed a simple analytical model for the computation of the CiS from 
textures. The model exploits the fact that in this scenario the microwave sky will show evidence 
of spots due to perturbations in the effective temperature of the photons resulting from the 
non-linear dynamics of concentrations of energy-gradients of the texture field. The model of 
course does not aim to replace the full range numerical simulations but just to show overall 
features predicted by textures in the CMB anisotropies. In fact the model leaves free a couple 
of parameters that are fed in from numerical simulations, like the number density of spots, u, 
the scaling size, d s , and the brightness factor of the particular spot, a k , telling us about its 
temperature relative to the mean sky temperature. 

Texture configurations giving rise to spots in the CMB are assumed to arise with a constant 
probability per Hubble volume and Hubble time. In an expanding universe one may compute 
the surface probability density of spots 

dP = N{y)dydQ, with N(y) = _ 8ul ^ ( 2 v/s - l) 2 , (1) 



1 The best place to learn about these missions and to follow the developments are, respectively, the sites 



http://map.gsfc.nasa.gov/ and http://astro.estec.esa.nl/SA-general/Projects/Cobras/cobras.html 



where Q stands for a solid angle on the two-sphere and the time variable y(t) = log 2 (t /i) 
measures how many times the Hubble radius has doubled since proper time t up to nowQ. 

In the present context the anisotropics arise from the superposition of the contribution 
coming from all the individual spots Sk produced from y\ s up to now, and so, AT/T = 
12k a kSk(0k,y), where the random variable stands for the brightness of the hot/cold k- 
th spot with characteristic values to be extracted from numerical simulations 6 ''. Sk(9k, y) is the 
characteristic shape of the spots produced at time y, where 9k is the angle in the sky mea- 
sured with respect to the center of the spot. A spot appearing at time y has typically a size 
6 s (y) ~ d s 6 hor (y), with 8 hor (y) the angular size of the horizon at y, and where it follows that 
9 s (y) = arcsin ( 2 »/3-i ) ■ Textures are essentially causal seeds and therefore the spots induced by 
their dynamics cannot exceed the size of the horizon at the time of formation, hence d s < 1. Fur- 
thermore the scaling hypothesis implies that the profiles satisfy Sk(0k,y) = S(8k/0 s (y)). From 
all this it follows a useful expression for the multipole coefficients, a™ = J2k a kSi(y)YJ n *(jk), 
with Sl(y) the Legendre transform of the spot profiles. 

At this point the C/s are easily calculated 5 -* . As we are mainly concerned with the three- 
point function we go on and compute the angular bispectrum predicted within this analytical 
model, which we find to be 

(W) = (« 3 ) / dyN{y)S\y)S\y)S\y) J d^Y^{^)Y^^)Y^ (7). (2) 

(a 3 ) is the mean cubic value of the spot brightness. 

Having the expression for the bispectrum we may just plug it in the formulae for the 
full mean three-point temperature correlation function 7 ). To make contact with experiments 
however we restrict ourselves to the collapsed case where two out of the three legs of the three- 
point function collapse and only one angle, say a, survives (this is in fact one of the cases 
analyzed for the four-year COBE-DMK data 8 -*). 

The collapsed three-point function thus calculated, (C 3 (a)), corresponds to the mean value 
expected in an ensemble of realizations. However, as we can observe just one particular real- 
ization, we have to take into account the spread of the distribution of the three-point function 
values when comparing a model prediction with the observational results. This is the well- 
known cosmic variance problem. We can estimate the range of expected values about the 
mean by the rms dispersion a^ v (a) = (Cf (a)) — (C^ia)) 2 . We will estimate the range for the 
amplitude of the three-point correlation function predicted by the model by (C 3 (a)) ±<7cv(a). 

It has been shown 6 ) that spots generated from random field configurations of concentrations 
of energy gradients lead to peak anisotropies 20 to 40 % smaller than those predicted by the 
spherically symmetric self-similar texture solution. These studies also suggest an asymmetry 
between maxima (a max ) and minima (a min ) of the peaks as being due to the fact that, for 
unwinding events, the minima are generated earlier in the evolution (photons climbing out of 
the collapsing texture) than the maxima (photons falling in the collapsing texture), and thus 
the field correlations are stronger for the maxima, which enhance the anisotropies. 

3 Results 

Let us now compute the predictions on the CMB non-Gaussian features derived from the 
present analytical texture model. One needs to have the distribution of the spot brightness 
{a k } in order to compute the mean values (a n ). It is enough for our present purposes to take for 
all hot spots the same > and for all the cold spots the same a c < 0. Then the (a n ) needed 

2 e.g., for a redshift z\ s ~ 1400 at last scattering we have y\ s ~ log 2 [(1400) 3 / 2 ] ~ 16. 
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Figure 1: The 'pseudo-collapsed' three-point function (in units of 10 4 /iK 3 ) as computed from 
the analysis of the four-year COBE-DMR data [data courtesy of G. Hinshaw and the COBE 
team]. The points plotted are the analogue of our ((73(a)). Error bars represent instrument 
noise while the grey band represents the rms range of fluctuations due to a superposition of 
instrument noise and cosmic variance. 

can readily be obtained in terms of (a 2 ) and x = (a)/(\a\). We fix (a 2 ) from the amplitude of 
the anisotropies according to four-year COBE-DMR 8 ). The other parameter, x, that measures 
the possible asymmetry between hot and cold spots, we leave as a free parameter. 

We consider the COB£HDMR window function and, in order to take into account the partial 
sky coverage due to the cut in the maps at Galactic latitudes |6| < 20°, we multiply acv by 
a factor ~ 1.5 in the numerical results (sample variance)^. Let us now compare with the 
data: Subtracting the dipole and for all reasonable values of the asymmetry parameter x, the 
data falls well within the (C 3 (a)) ± acv( a ) band, and thus there is good agreement with the 
observations. However, the band for Gaussian distributed fluctuations (e.g., as predicted by 
inflation) also encompasses the data well enough, and it is in turn included inside the texture 
predicted band. The fact that the range of expected values for the three-point correlation 
function predicted by inflation is included into that predicted by textures for all the angles, 
and that the data points fall within them, makes it impossible to draw conclusions favoring one 
of the models. 

It is well known that the largest contribution to the cosmic variance comes from the small 
values of I. Thus, no doubt the situation may improve if one subtracts the lower order multipoles 
contribution, as in a i m in = 10 analysis 9 ^. In Figure 1 we show the analysis of the four-year 
COBE-DMR data evaluated from the 53 + 90 GHz combined map, containing power from the 
t = 10 moment and up. It is apparent that the fluctuations about zero correlation (i.e., no 

3 In the analysis of the data, the method used for computing the uncertainty is to generate 2000 random 
skies with HZ signal + noise, then to compute the three-point function of each realization on the cut sky after 
subtracting a best-fit multipole. Hence, the result automatically includes sample variance. 




Figure 2: Expected range of values for the three-point correlation function for the texture 
model with asymmetry parameter x = 0.4 (grey band). Also shown is the expected range for 
inflationary models (black band). Both bands include the ~ 50% increase in acv due to the 
sample variance. All multipoles up to I = 9 have been subtracted. The right panel shows a 
zoomed fraction of the same plot. 

signal) are too large for the instrument noise to be the only responsible. These are however 
consistent with the range of fluctuations expected from a Gaussian process (grey band). 

What we want to see now is whether our analytical texture model for the three-point 
function 10 ^ can do better when compared with the data. Figure 2 shows the collapsed three- 
point function (C 3 (a)) (solid line) and the grey band indicates the rms range of fluctuations 
expected from the cosmic variance. Also shown is a black band with the expected range for 
inflationary models (no instrument noise included). The bands do not superpose each other 
for some ranges of values of the separation angle for the value of x — 0.4 considered, what 
means that measurements in that ranges can distinguish among the models. The value of the 
parameter x considered is quite in excess of that suggested by simulations 6 ^ , but was chosen to 
show an example with a noticeable effect. 

In Figure 3 we show the result of combining the previous figures, confronting the actual data 
with the curves predicted by the texture model. From these figures one may see qualitatively 
by eye that (for some ranges of the angular separation better than for others, of course) the 
data seems to follow 'approximately' the trend of the texture curves. While many of the data 
points fell outside the Gaussian band (Figure 1), most of them are now inside the grey band 
in Figure 3 (left panel). Moreover, while we vary the x values from 0.4 down to 0.1 (the actual 
value suggested by texture simulations) we see that more and more points (with error bars) get 
inside (or touch) the grey band. Can this be just by chance? Or is it there something worth of 
further study? In order to answer these questions one ought to quantify more the analysis by 
using a \ 2 statistics for the model and data, and comparing it to the Gaussian case, e.g. 11 ) 

X 2 = £(A* " (Cs) a )(M-%,(D p - (C s ) fi ), (3) 

a/3 

with D a the COBE-DMR three-point function and M is the covariance matrix of the analytical 
model. It might be that the data picks out a preferred value for the asymmetry parameter x. 
Work in this direction is currently under way. 




Figure 3: Combined four-year COBE-DMR data and collapsed three-point correlation function 
predicted by the analytical texture model (as in previous figures). Left panel: for a somewhat 
exaggerated value of the asymmetry parameter x = 0.4. Right panel: for the value x ~ 0.1 
suggested by texture simulations. 
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